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ABSTRACT 

The  interaction  between  the  mechanical  deformation  and  the 
thermal  state  of  solids  is  studied.  Principles  of  classical 
thermodynamics  are  used  to  calculate  temperature  changes  arising 
from  the  isentropic  deformation  of  elastic  solids.  The  theory  of 
irreversible  thermodynamics  is  used  to  predict  temperature  changes 
due  to  the  adiabatic  deformation  of  viscoelastic  solids. 

It  is  found  that  such  temperature  changes  are  small  and  that 

the  effect  of  thermomechanical  coupling  may  usually  be  neglected 

I * 

in  the  non-oscillatory  deformation  of  linear  elastic,  finite 
(rubber)  elastic,  and  linear  viscoelastic  solids.  Data  are  pre- 
sented for  several  metals,  rubbers,  and  plastics. 


» • 


INTRODUCTION 

In  the  analysis  of  mechanical  deformation,  thermal  effects  are 
usually  not  considered  beyond  the  occasional  noting  of  the  temperature 
of  the  material  at  the  onset  of  testing.  Because  the  mechanical  prop- 
erties of  the  material  can  be  essentially  independent  of  the  temperature 
over  large  temperature  ranges,  this  approximation  is  usually  justifiable. 
Even  when  the  material  properties  exhibit  significant  thermal  depen- 
dence, the  deformation  process  can  cause  such  a sma^l  change  in  tempera- 
ture that  its  effects  may  be  ignored. 

There  are  some  cases,  however,  especially  in  the  realm  of  visco- 
elastic materials,  in  which  the  temperature  change  and  the  ensuing 
variation  of  mechanical  properties  can  be  important. 

We  shall  concern  ourselves  with  essentially  static  loadings,  such 
as  simple  extension,  shear,  or  bulk  compression  of  solids.  These 
"static"  processes  will  be  referred  to  as  transient.  Changes  in  the 
temperature  of  a body  which  result  from  its  deformation  under  such 
loadings  will  be  calculated.  Knowledge  of  the  effect  of  temperature 
upon  the  mechanical  properties  then  enables  us  to  gauge  the  importance 
of  thermomechanical  coupling  in  a given  problem. 

The  following  chapters  will  deal  with  the  transient  loadings  of 
linear  elastic,  rubber  elastic,  and  linear  viscoelastic  solids.  We 
will  not  consider  here  dynamic  processes  such  as  vibration  or  shock 
loading. 


A 
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CHAPTER  ONE 


LINEAR  ELASTIC  DEFORMATIONS 


The  effect  of  thermomechanical  coupling  in  linear  elastic  solids 
is  easily  derived  from  the  basic  equations  of  thermodynamics  and  thermo- 
elasticity. Results  in  this  field  are  already  well  established  (see 
for  example  Nowacki  [1]  and  Biot  [2]), and  the  following  derivation  is 
presented  primarily  as  a basis  and  as  a source  of  comparison  for  subse- 
quent developments. 

In  the  ensuing  equations,  - uantities  such  as  internal  energy,  U , 
heat,  Q , work,  W , entropy,  S , et  cetera,  are  measured  per  unit  of 
original  volume  of  an  infinitesimal  body.  We  write  the  first  law  of 
thermodynamics  as^ 


dU  = dQ  + dW 


The  absolute  temperature  of  the  body  is  denoted  by  T . From  the  second 
law  of  thermodynamics,  applied  to  reversible  processes,  we  write: 


dQ  = TdS 


The  work  done  on  the  body,  per  unit  of  original  volume,  is  expressed  by^ 


dW  = a . .de . 
ij  ij 


1 dW  is  the  total  increment  of  work  done  on  (not  by)  the  body,  and  it 
includes  work  done  by  atmospheric  pressure  and  any  other  pre-loads. 
These  are  departures  from  usual  thermodynamic  convention.  dQ  denotes 
the  increment  of  heat  added  to  the  body. 

2 Summation  convention  will  be  used,  implying  summation  over  repeated 

free  indices.  Thus  3 3 

o .de  = l Z a.  .de.., 
ij  i-1  j-1  ^ « 


2 


where  o„  and  e„  are  the  conventional  stresses  and  infinitesimal 

3 

strains  of  classical  elasticity. 

Using  the  definition  of  the  Helmholtz  free  energy, 

A = U - TS  , (4) 

entropy  may  be  expressed  as  a function  of  temperature  and  the  strain 
tensor,  S = S(T,e)  . From  eqs.  (1-4),  we  have 

dA  = a. .dc  . - SdT  (5) 

ij  ij 

4 

Thus : 


3A 

de.  . 
ij 


a . . 
ij 


(6a, b) 


It  follows  from  the  commutivity  of  partial  differentiation  and  the 
regular  behavior  of  A that 


i'A 

cEj 


da.  . 

hi 

8T 


'T,e' 


From  eq.  (2)  and  the  definition  of  heat  capacity,  we  write 


(7) 


3 To  be  precise,  these  can  be  interpreted  as  the  second  Piola-Kirchof f 
stress  (force  per  unit  of  undeformed  area)  and  the  Lagrange  strain 
(attributed  to  Cauchy.)  In  the  infinitesimal  case,  these  measures 
differ  insignificantly  from  various  other  common  measures.  For 
detailed  discussion  of  these  points,  see  Truesdell  and  Noll  [3]  and 
Rivlin  [4]. 

4 The  notation  | „ implies  that  all  elements  of  £ are  held  fixed  with 
the  exception  — of  the  one  indicated  in  the  partial  differential 
operator. 


3 


8S_ 

3T 


(8) 


where  pQ  is  the  material  density  in  the  undeformed  state  and  Cy 
the  heat  capacity  per  unit  mass  at  constant  volume  or,  to  be  more 
precise,  constant  strain. 

The  partial  derivative  on  the  right-hand  side  of  eq.  (7)  may  be 
evaluated  directly  from  consideration  of  the  therraoelastic  stress- 
strain  laws.  Letting  y ana  X denote  the  isothermal  Lam£  constants, 
a the  volumetric  coefficient  of  thermal  expansion,  and  8 the  increase 
in  temperature  from  the  strain-free  state,  strain  is  given  by 


- —f  > A 

ij  2ul°ij'°o  ' 2u(2B+3») 


<V  - %, , > «1,  + 3esn  • <»> 


where  a is  the  pre-loading,  the  loads  existing  in  the  strain-free 

°ij 

state.  The  most  common  example  of  pre-loading  is  atmospheric  pres- 

sure. Inverting  these  relations,  we  have 


°ij  = V.  + 2u£ij  + X£kk6ij  - aK06ij  (10) 

v;here  K = X + — j-  is  the  isothermal  bulk  modulus. 

Writing  entropy  in  terms  of  its  partial  derivatives,  and  substi- 
tuting from  eqs.  (7),  (8),  and  (10): 


dS  = dT  + aK6  .de,  . 

T ij  ij 


Thermomechanical  coupling  effects  are  usually  ignored  in  non- 


oscillatory  situations.  An  argument  often  cited  in  support  of  this 


4 


approximation  is  one  which  we  shall  refer  to  here  as  the  "equivalent 

heat"  argument.  Stated  simply,  it  states  that  if  one  were  to  supply 

the  thermodynamic  system  (the  infinitesimal  body)  with  an  amount  of 

heat  equivalent  to  the  mechanical  work  done  on  the  system,  the  resulting 

change  in  temperature  of  the  system  would  be  an  upper  bound  for  the 

actual  change  in  temperature.  (The  change  in  temperature  predicted  by 

this  argument  is  denoted  by  dT  . ) Thus: 

eqh 


dW  = dQ  = p C.  dT  . 

o V eqh 


and  from  eq.  (3) , 


o . .de . , 
r = _JJ U 

eqh  PC 


dTeqh  is  compared  with  the  result  of  the  thermodynamic  analysis, 
denoted  by  ^ * For  an  adiabatic,  reversible  process,  dS  = 0 
and  eq.  (11)  becomes: 


- ctTK  , 

poCV  lj  ^ 


- aTK 

P C d£kk 
o V 


Integrating  eqs.  (13)  and  (14),  making  the  approximations  that 
a and  K are  constant,  is  linear  in  temperature,  and  AT  is 
small  with  respect  to  T , we  find  that: 


- aTK 


r 


Comparison  of  eqs.  (15)  and  (16)  shows  that  AT  , is  a poor  measure 

eqn 

of  AT  , the  actual  temperature  change.  AT  , is  dependent  solely 
thm  thm 

upon  the  volume  change  and  is  thus  identically  zero  for  any  equivo- 
luminal  deformation.  ATg  ^ > however,  is  proportional  to  the  incre- 
ment in  strain  energy  and  is  affected  by  shear  deformations.  In  such 
cases,  ATg  ^ does  indeed  serve  as  an  upper  bound  for  the  change  in 
temperature,  but  it  is  misleading  to  place  an  upper  bound  on  a quantity 
which  must  be  zero. 

In  the  case  of  hydrostatic  compression,  AT  ^ can  fall  far 
short  of  AT  . As  an  example,  consider  the  compression  of  a speci- 
men by  a pressure  of  n atmospheres.  Denoting  one  atmosphere  of 

pressure  by  P , then  a = - P6..  , Ao..  = - nP6 . . , and 

Oy  iJ  iJ  ij 

Ae, . = - nP6../3K  . Substitution  into  eqs.  (15)  and  (16)  leads  to: 
ij  ij 


AT 

AT 


thm 

eqh 


2aTK 
(2+n) P 


(17) 


The  ratio  can  be  quite  large,  as  is  seen  in  the  following  tabulation: 


TABLE  1 : 3 Sample  Values  of  (aTK/P) 


p = 14.7  psi  T = 

298°K 

Material 

a,(°K_1) 

K, (psi) 

aTK/P,  (— ) 

mild  steel 

3.6xl0~5 

24.6xl06 

1.8xl04 

polyisoprene (ebonite  rubber) 

1.9xl0~4 

6.0xl05 

2.3xl03 

polyethylene (low  density) 

7.4xlO-4 

4.8xl05 

7. 2xl03 

5 From  [5],  table  1-59;  [6],  p.  708 ; [ 7 ] , sec  VI;  [8],  table  9.1  . 


Even  in  the  case  of  ebonite,  the  most  favorable  case  presented,  the 

ratio  of  AT  to  AT  does  not  reach  unity  until  the  material 
thm  eqh 

is  compressed  by  almost  4600  atmospheres.  Pearson  [9]  presents  a 

calculation  for  a cube  of  steel  compressed  from  an  initially  unstressed 

state.  In  his  example,  the  ratio  of  AT  to  AT  is  unbounded  as 

thm  eqh 

the  applied  pressure  approaches  zero,  as  may  be  seen  by  examination 

of  eqs.  (15)  and  (16)  with  o » 0 . 

°i) 

Despite  the  failure  of  ^ to  predict,  or  even  to  provide  an 

upper  bound  on,  the  thermal  effects  of  deformation,  these  effects  are 
still  negligible  for  infinitesimal  elastic  deformation,  as  is  later 
seen  in  Chapter  Four. 
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CHAPTER  TWO 


FINITE  ELASTIC  DEFORMATIONS 


Since  many  polymers  exhibit  rubber  elasticity,  we  perform  a similar 
analysis  for  finite  elasticity.  We  will  deal  only  with  principal  stresses, 
and  will  use  extension  ratios  and  true  stress  (force  per  unit  of  deformed 
area,  sometimes  called  Cauchy  stress)  defined  as  follows:8^ 

Xi  = Li/Lo(i)  = 1 + £i(i)  (18) 


ti  fi/LjLk  °i(i)/AjXk  * 


In  terms  of  these  variables,  the  first  law  of  thermodynamics 
becomes,  for  reversible  processes: 


dU  = TdS  + (V/V  ) t . dLnX . 

o 1 1 


Eq.  (5)  becomes: 


dA  = (V/V  ) t dLnA  - SdT 
o i i 


We  proceed  as  in  the  infinitesimal  case,  noting  that  V = ^X^^X^  ’ 

(T\  8 

X = (A^,  X^,  X^)  J implies  constant  volume. 


6 The  extension  ratio  is  closely  related  to  the  natural,  or  logarithmic, 

strain;  = LnX^  . 

7 Numerical  indices  and  indices  hidden  in  parentheses  are  not  subject 
to  implied  summation. 


8  T used  as  a superscript  to  a vector  or  a matrix  will  denote  the 
transpose  of  that  quantity.  T used  as  a scalar  will  continue  to 
denote  absolute  temperature. 


8 


As  in  eq  . (8)  , 


3S 

3T 


' “oV1 


(22) 


Similarly  to  eqs.  (6), 


3A 

3T 


= - S 


3A 
3A . 


T,r 


V(i) 


(23a, b) 


Thus : 


3S 
3A . 


Vd)  3T 


(24) 


From  eqs.  (22)  and  (24)  we  have 


(25) 


(26) 


(27) 


(28) 


AT 


thm  p V 
o 


o f V 3ti 

oCV  J X(i)  3T 


dA. 


AT 


eqh 


' '.V>  / xa>  tldil 


(29) 


At  this  point  we  consider  possible  forms  of  the  function 
Ci  = In  most  references,  it  is  assumed  that  stress  varies 

as  a linear  function  of  temperature  in  strain  control  experiments. 

This  assumption  is  well  verified  by  experimental  results  (Meyer  and 
Ferri  [10]  and  Anthony,  Caston,  and  Guth  [11]  as  cited  in  Treloar  [12]). 
As  a further  motivation,  we  shall  show  that  this  linear  dependence  is 
a direct  result  of  another  physical  observation,  that  the  heat  capac- 
ity at  constant  volume  is  independent  of  the  deformational  state  over 

large  ranges  of  deformation,  that  is  to  say, 


CV(A,T)  = Cy(T) 


(30) 


Taking  eqs.  (22)  and  (24)  and  invoking  the  commutivity  of  partial 
differentiation,  derivation  of  the  second  partial  derivate  of  S with 
respect  to  X and  T requires  that 


P 

o 

T 


-V 


32t, 


T»X  " 


V A,,.  __2 
o (i)  3T 


(31) 


Eq.  (31)  and  the  form  of  eq.  (30)  require  that 


(32) 


which,  integrated  twice,  becomes: 


10 


ti(A,T)  = hi(A)T  + g.(x)  (33) 

Conversely,  assumption  of  eq.  (33)  leads  to  eq.  (30).  Using  the  form 
of  eq.  (33),  eqs.  (28)  and  (29)  become: 


v c„  / ; 

° o VJ 


thm  p 


AT  . = AT  L 
eqh  thm 


■p-  h (A)dA 

A (i)  1 ' 1 


+ TvT  / T~"  8i(x)  dXi 
Po  o V J X(i)  1 - 1 


As  a sample  calculation  using  eqs.  (34)  and  (35),  we  consider  the 
uniaxial  extension  of  natural  rubber  and  use  data  given  in  Treloar 
[12].  Using  the  simple  extension  law  for  a Mooney  material  and  the 
form  of  Eq.  (33)  as  guides,  the  data  presented  may  be  reasonably  well 
fitted  by: 


t ( A ,T)  = [ . 01 (A2  - A-1)-. 001]  T + 1. 7 (A-A  2)  +0.4  (Kg/cm2)  (36) 


Letting  V = and  T = 300°K  , tables  of  P^^T 

p C (AT  - AT  ) are  presented, 

o V eqh  thm 


TABLE  2a:  p CTI  AT 

o V thm 


TABLE  2b:  p C„(AT  L - AT.,  ) 

o V eqh  thm 


. Oui 

.01 

.1 

1. 

.01 

.1 

1. 

* 

* 

.014 

2.8 

1.0  * 

* 

.061 

1.3 

.00055 

.0058 

.092 

3.2 

1.1  .00078 

.0080 

.093 

1.4 

.0013 

.013 

.16 

3.9 

1.2  .00105 

.0106 

.115 

1.5 

.0030 

.030 

.32 

5.0 

1.5  .0015 

.015 

.15 

1.7 

,0051 

.051 

.53 

6.9 

2.0  .0017 

.017 

.17 

1.7 

,0086 

.087 

.87 

10.2 

3.0  .0018 

.018 

00 

»H 

1.8 

(table  values  in  Kg/cm2) 
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No  calculations  are  made  at  points  marked  with  (*)  because  the  results 
would  be  far  too  sensitive  to  the  constants  chosen  in  eq.  (36)  to  be 
meaningful. 

In  the  case  of  uniaxial  extension,  ^ consistently  overesti- 

mates the  temperature  rise.  This  overestimation  varies  from  around  20% 
to  over  400%  in  the  examples  calculated,  and  is  seen  to  be  more  reason- 
able for  higher  initial  and  subsequent  extensions.  Again,  the  actual 
temperature  rise  is  seen  to  be  negligible.  For  the  most  extreme  case 
shown,  extension  of  a specimen  from  X = 3 to  X = 4 , thermodynamics 
predicts  an  increase  of  0.65°C  (AT^^  is  0.76°C.)  In  rubber  elasticity, 
such  effects  may  be  ignored. 


12 


CHAPTER  THREE 

LINEAR  VISCOELASTIC  DEFORMATIONS 

The  thermodynamic  analysis  of  viscoelastic  materials  is  far  more 
complicated  than  that  of  elastic  media.  In  the  viscoelastic  case, 
stress  and  strain  are  no  longer  simple  functions  of  each  other,  but 
can  instead  be  related  through  differential  operators  or  through  more 
general  functionals.  Furthermore,  deformation  processes  are  irrevers- 
ible. Comprehensive  discussion  of  these  problems  may  be  found  in 
several  texts,  and  reference  is  made  here  to  Fliigge  [13],  Lee  [14],  and 
Christensen  [15]. 

In  this  work,  we  shall  deal  with  stress-strain  relations  in  the  form 
of  hereditary  integrals  over  the  entire  previous  history  of  the  material. 

We  shall  restrict  ourselves  to  consideration  of  viscoelastic  materials 
which  can  be  described  as  linear  in  the  sense  that  a Boltzmann  superposi- 
tion may  be  employed  for  complicated  or  continuously  varying  stress  or 
strain  histories.  Discussion  of  this  technique  may  be  found  in  references 
[13-15],  and  further  comments,  including  some  examples  of  limitations  and 
exceptions  to  this  linearity,  are  available  in  Staverman  and  Schwarzl  [16], 
Leaderman  [17],  and  Matsuoka,  Alolsio,  and  Blair  [18]. 

Further  complication  is  introduced  by  the  fact  that  viscoelastic 
material  properties  exhibit  marked  temperature  dependence.  Polymer  chem- 
ists have  long  noted  this  behavior  and  have  searched  for  good  macroscopic 
descriptions  as  well  as  microscopic  explanations  for  the  phenomena. 

Ferry  [19]  and  Tobolsky,  in  conjunction  with  several  co-workers  [20-23], 
provide  many  data  showing  a particular  relationship  between  the  tempera- 
ture and  the  time  history  of  mechanical  deformations.  Schwarzl  and 


13 


Staverman  [24]  discuss  the  physics  of  this  relationship  and  provide  the 
term  "thermorheologically  simple"  to  denote  the  class  of  materials 
exhibiting  it.  Morland  and  Lee  [25]  have  generalized  the  time-temperature 
relationship  to  include  the  effects  of  continuously  varying  temperature 
history.  Many  polymers  can  be  modeled  as  thermorheologically  simple. 

This  work,  however,  deals  with  small  temperature  excursions  from  the 
equilibrium  state,  and  the  question  of  thermorheological  simplicity 
need  not  be  considered. 

Rivlin  [26]  and  Green  and  Rivlin  [27]  use  functionals  to  express 
viscoelastic  mechanical  properties  and  they  examine  the  form  which  these 
functionals  may  assume.  Coleman  [28]  applies  this  procedure  to  the 
remaining  thermodynamic  variables.  The  same  approach  is  developed  in  a 
more  phenomenological  manner  by  Hunter  [29],  Pipkin  [30],  Breuer  and 
Onat  [31],  Christensen  and  Naghdi  [32],  and  Cost  [33].  In  their  work, 
material  laws  are  expressed  in  single  and  multiple  integral  form. 

The  dissipation  inherent  in  viscoelasticity  necessitates  the  consid- 
eration of  non-equilibrium,  irreversible  thermodynamics.  Two  of  the 

classical  texts  in  this  field,  providing  reference  to  the  pioneering 
efforts  and  discussion  of  the  basic  principles  involved,  are  those  of 
De  Groot  [34]  and  Prigogine  [35]. 

In  the  application  of  irreversible  thermodynamics  to  viscoelasticity, 
Biot  [36,37]  expresses  entropy  in  "near-equilibrium"  states  as  a quad- 
ratic expansion  in  an  indefinitely  large  number  of  "hidden  coordinates" 
and  proceeds  to  derive  stress-strain  laws  analogous  to  spring-dashpot 
models  of  viscoelasticity.  Eringen  [38]  and  Schapery  [39]  apply  this 
procedure  to  several  problems  in  continuum  mechanics. 
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In  1966,  Valanis  [40]  published  a paper  showing  that  solutions  to 
the  equations  resulting  from  the  hidden  variable  approach  may  be  expressed 
as  integral  operators  of  the  form  employed  in  the  phenomenological  approach. 


! 


In  this  chapter,  thermodynamical  equations  will  be  formulated  in  terms 
of  hidden  coordinates,  ;s  in  reference  [39],  and  will  be  solved  in  the  man- 
ner suggested  by  Valanis  [40].  The  solutions  will  then  be  used  to  prove 
an  important  symmetry  property;  this,  in  turn,  will  enable  us  to  apply  the 
results  to  the  question  of  thermomechanical  coupling  effects. 

Keeping  our  previous  notation  for  entropy,  internal  energy  (both  per 

unit  original  volume X and  temperature,  we  now  denote  stress  and  strain  by 

six-element  column  vectors,  a and  q^  , respectively,  such  that  dW  = 

T 

a dq^  . A column  vector  of  "hidden  variables,"  q 2 , is  defined  as  the 
vector  whose  elements  consist  of  all  the  variables  necessary  to  complete 
the  definition  of  the  thermodynamic  state.  The  unknown  number  of  elements 
of  g will  be  denoted  by  N . As  before,  our  thermodynamic  system  is 
an  infinitesimal,  homogeneous  body  having  uniform  temperature. 

In  the  manner  of  De  Groot  [34],  we  define  the  quantity  dSj  , the 
internal  entropy  production,  by: 


dSj  = dS  - dQ/T  (37) 

We  express  the  entropy  as  a function  of  the  state  variables,  U , q^  , 
and  q^  , and  introduce  the  convenient  notation: 


S 


No  work  can  be  done  if  the  deformations,  , are  held  fixed,  and 

there  are  no  internal  irreversibilities  associated  with  processes  for 
which  dq^  = 0 = dq2  (see  Schapery  [39]).  Therefore,  from  eqs.  (1) 
and  (2): 


Dividing  both  sides  of  eq.  (45)  by  a small  time  increment,  we  write: 


• T T 

TSi  = ?i  Si  + h 32 


where  a superposed  dot  denotes  differentiation  with  respect  to  time 
or,  as  later  used  under  an  integral,  with  respect  to  the  time-like 


argument. 


We  now  seek  more  information  about  the  function  S^.  . The 
T T 

coefficients,  and  J^  , in  eq . (46)  can  depend  upon  the  state 

variables,  T , q^  , and  q^  , and,  additionally,  upon  the  applied 
stresses,  a . If  we  postulate  a relationship  between  these  quantities 
and  the  rates  of  change  of  the  system,  q , q , and  0,  we  can  write 
Sj  as  a function  of  these  latter  variables  directly.  Accordingly,  we 
shall  expand  S in  a truncated  Taylor  series  about  an  equilibrium 
state  (in  which  0 = 0 , q^  = 0 = q ) , neglecting  all  terms  of  higher 
than  second  order. 


SI  = “o  + *l\  + *2\  + V + + W + 


+ I?2T^82T<?1  + e3?2^  + «?l\  + ^2^2^  + 2^3®2 


* and  ^3  are  constant  matrices,  0 and  B^  being 
symmetric.  This  expansion  may  be  simplified  by  noting  that,  accord- 


ing to  the  second  Jaw  of  thermodynamics,  = 0 in  an  equilibrium 

state;  thus  a = 0 . The  second  law  also  requires  that  S_  never 
o I 

be  negative.  Since  the  q^'s  anc*  9 may  be  arbitrarily  small  and  are 


TT*—  — 


not  restricted  as  to  sign,  this  requires  that  a,  = 0 = a„  , that 

~1  ~ ~ 2 

s 0 , and  that  0^  and  0^  be  positive  semi-definite.  Since 
the  elements  of  are  all  associated  with  internal  dissipation, 

0^  must  actually  be  positive  definite.  Eq.  (47)  now  reduces  to 


" i[?lT(Bl?l  + B2?2  + !lB)  + ?2T<62T5l  + S3?2  + I2B)  + 


+ (-lT^l  + ll'q2  + 


(48) 


Comparison  with  eq.  (46)  reveals  that  the  three  sets  of  terms  in 

T 

parentheses  are  the  three  "irreversible  forces"  associated  with  q^  » 

. T 

^2  , and  0 . We  note,  however,  that  the  third  term  is  zero,  since 

0 does  not  appear  explicitly  in  eq.  (46).  This  is  reasonable  since  an 
increase  in  temperature  of  the  system,  if  not  accompanied  by  any  changes 
in  the  mechanical  or  "hidden"  coordinates,  does  not  cause  internal  dissi- 
pation. This  further  reduces  eq.  (48)  by  requiring  that  = 0 and 

that  Yj  = 0 = Y2  • We  now  write  the  "irreversible  forces"  as: 

h = K(6lSl  + (49a) 

*2  = iTo(62TSl  + e3S2}  (49b) 


where  T becomes  Tq  in  keeping  with  the  truncation  of  the  power 
series  (eq.  (47))  to  second  order,  and  solutions  for  the  involving 
antisymmetric  matrices  are  suppressed  in  accordance  with  the  symmetry 
requirements  of  Onsager's  principle  [41,42].  It  should  also  be  noted 
that  the  quadratic  terms  in  eq.  (47),  resulting  ij.  n our  assumption  of 
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a relationship  between  the  q^'s  and  9 and  the  state  variables, 
reflect  Onsager's  principle  of  linear  relationship  between  "forces 
and  fluxes."  Onsager's  principle  is  restricted  to  near-equilibrium 
states;  our  derivation  is  also,  due  to  the  truncation  of  the  Taylor 
series  after  the  quadratic  terms. 

We  now  introduce  the  free  enthalpy  (Gibbs'  function): 

G = U - TS  - goTq  (50) 

where  g^  is  the  stress  in  the  unstrained  state.  Differentiating 
eq.  (50)  and  substituting  from  eqs.  (41-43): 

dG  = - SdT  - [ao  + T(13S)]Tdq1  - T(2aS)Tdq2  (51) 

so  that 


3G 


- - [oo  + T(  9S) ] 


T’S£’S2 


3G 

3?2 


= - t(  as) 


aG 

ai 


T’?l,q2 


= - S 


qrq2 


(52a) 


(52b, c) 


If  the  free  enthalpy  is  expanded  about  an  unstrained  equilibrium 
state  (in  which  0 = 0 , q^  = 0 = q^) , we  obtain: 


G = Go  + a30  + a^  + a2Tq2  - f^e2  + ^/(B^  + 8^)  + 
+ 2?2T(B2Tqi  + B3q2}  + (^lTql  + ^2Tq2)e 


(53) 


J 


f 


j 
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where  and  B^  are  symmetric.  It  may  be  shown  (see  ref.  [39]) 

that  G _>  Gq  for  any  nearby  states  in  which  0=0.  Therefore, 

a^  = 0 = a^  and  B^  and  B^  must  be  positive  semi-definite. 

Irreversibility  further  requires  that  B^  be  positive  definite. 

2 

The  negative  sign  is  assigned  to  the  0 term  for  later  convenience. 
Differentiating  eq.  (53), 


9^  ' - T<  - So  ' Bl?l  + B2?2  + h6 


(54a) 


I?  - - t(2ss)  - b2t?1  + b3,2  + b2e 


9?2 


(54b) 


It  = ~S  = a3  + ^1  + ' b36 


(54c) 


Hence,  from  eqs.  (44),  (49),  and  (54), 


2 - % ■ iTo<Bl?l  + B2?2>  + BlSl  + B2?2  + h6 


(55a) 


o ■ I V62T5l  + B3?2>  + B2T8l  + B3?2  + h® 


(55b) 


We  note  that  B^  and  0^  may  be  simultaneously  diagonalized 
by  a non-singular  matrix,  V (see  Franklin  [43]).  Since  both  matrices 
are  positive  definite,  we  may  further  require  that  V be  such  that 


VTB3V  = A and  V^V  = I 


(56a ,b) 
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where  I is  the  N x n identity  matrix  and  A is  a diagonal 


matrix  of  positive  constants.  If  we  define  the  transformation  of 


Analogous  to  a result  in  Valanis  [40] 


Integration  by  parts  yields  a more  desirable 


form  of  this  equation 


where 


and  E commutes  with  A and  A 


The  form  of  eq.  (60)  is  of  great  significance.  We  have  started 
with  N "hidden  coordinates",  , where  N is  an  unknown,  arbi- 

trarily large  number.  After  rearranging  these  coordinates  in  a con- 
venient manner  (see  eq.  (57)),  we  find  that  they  are  simply  various 
complicated  expressions  of  the  previous  deformation  and  temperature 
histories.  Thus  the  current  state  of  the  thermodynamic  system  is 
prescribed  as  a functional  of  the  temperature  and  the  physical  defor- 
mation from  the  onset  of  activity  to  the  current  time.  Furthermore, 
the  fact  that  the  elements  of  A are  positive  shows  that  this  func- 
tional is  similar  to  the  relaxation  modulus,  familiar  in  isothermal 
viscoelasticity. 

Applying  eqs.  (57),  (58),  and  (60)  to  eq.  (55a),  we  find  that 


a - go  = [B][  - B2VA  V^lq  + [t^  - B^A'Vb^B  + 


+ K[B1  - 62WVl?l  + 


j [b2va-1 


B2V]Ex(t-x)V  b.,0dT  + 


- j [B2VA_1  - B2V]E  (t-T)[AVTB2T  - VTB2T]q1dT 


Referring  to  eq.  (54c)  and  taking  a^  = - S = - Sq 


application  of  eqs.  (57)  and  (60)  yields: 


22 


(63) 


s ~ s0  = tb2  VA  ViB2i  - b^]^  + [b3  + b2  VA  V b2]0  + 


T 

+ b2  V 


/ Ex(t-T 


) [VT32T  - A_1VTB2T]q1dT  + 


- '2 


rv  J Ex(t-T)A-1vTb2edT 


Eqs.  (62)  and  (63)  can  be  compared  with  those  resulting  from  the 
phenomenological  approach,  in  this  case  eqs.  3.9  and  3.10  in 
Christensen  [15].  Modifying  his  notation,  to  bring  it  in  line  with 
this  development,  we  write: 


a - a = 
' o 


S - S = 
o 


rt  rc 

I G(t-x,0)  q^dx  -/  <J>(0,  t— t)  9dr 
o *o 

/ <fiT(t-T,0)  q^dx  + / m(t-T,O)0 

o *0 


(64a) 


0dt 


(64b) 


Taking  the  q^  and  0 terms  in  eqs.  (62)  and  (63)  under  the  integral: 
we  make  the  association: 


-1„T„  T 


G(t-r)  = B - B VA  V B0  - [B.V  - B VA]  E (t-x)  . 
It  t 2 2 x 

• [VT62T  - A"Vb2T]  + - 62WT6^]6(T-t) 


(65a) 


t„a-i„t, 


m(t-t)  - b + b AVA  AVAb  - b XVE  (t-x)A-1VTb 


X 


(65b) 


4>(0,t-x)  - - b.  + B„VA  VAb,  - B VE  (t-x)A  V V + 
'1  2 -l  2 x '2 


(65c) 


+ B,VE  (t-t)V1b, 


T T T —1  T T T 

<t>  (t-T,0)  = - b + b2AVA  V1B21  + b2  VEx(t-x) 


[VTB2T  - A 1VTB2T1 


It  can  now  be  seen  from  eqs.  (65c&d)  that  <j>(0,t)  = <t>(t,0)  , 
a symmetry  property  which  is  assumed  in  many  works  [15,29,32,33] 
but  apparently  has  not  been  proven  previously. 

Another  interesting  feature  of  the  stress-strain  law,  eq.  (62), 
is  the  appearance  of  the  q^  term  outside  the  integrals.  This  term 
represents  the  dependence  of  the  stress  state  upon  the  instantaneous 
strain  rate,  which  is  now  seen  to  be  a natural  result  of  the  second- 
order  approximations  of  entropy  production  and  free  enthalpy.  The 
term  may  be  incorporated  under  the  strain-rate  history  integral  by 
use  of  the  Dirac  delta  function.  This  has  been  done  in  eq.  (65a)  to 
facilitate  comparison  of  eqs.  (62)  and  (64a).  It  would  seem  advanta- 
geous, however,  to  consider  this  term  separately  from  the  relaxation 
modulus.  Its  presence  precludes  any  instantaneous  elastic  response. 
Since  solid  materials  generally  exhibit  a measurable  instantaneous 
response,  the  coefficient  of  q^  is  customarily  taken  to  be  zero. 


Much  of  the  derivation  up  to  this  point  has  been  aimed  at 
deducing  this  symmetry  of  <|>  in  its  two  arguments.  With  this  accom- 
plished, and  with  the  machinery  at  hand,  we  address  the  question  of 
therraomechanical  coupling. 

Since  most  polymeric  materials  have  poor  thermal  conductivity 
and  since  small  temperature  changes  are  anticipated,  it  is  reasonable 
to  model  most  transient  deformation  tests  as  adiabatic  processes. 
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This  approximation  reduces  eq.  (37)  to  dS^.  = dS  , which  we  shall 
write  as 


ST  - S (66) 


From  eqs.  (48)  and  (54c),  we  then  have: 


2^?1  3lSl  + 2?1  e2?2  + ^2  B3?2^  “ b39  “ -1  ?1  ~ -2^2  (67^ 


After  application  of  eqs.  (56b)  and  (57), 


0 = 2?lTei?l  + ?lTe2V^  + 2^3  ~ b30  + blT^l  + b2^V^ 


(68) 


After  elimination  of  q , using  eq.  (58),  and  much  cancellation  of 
terms : 


• 2 T T T T T T 

0 = - b3  0 - - b2  W b2  9 + [b^  - b2iWi82i]q1  + 

o 


2TTT  4 TTT  2TT2 

- Tq-2  W B2  ?1  + "2  h VV  »2  V + “2^2  VV  V + 

o o 


(69) 


2T  TT  IT  TT  o t 

+ -2  Sl  B2  W B2  ?1  + [61  - 82W182T]q1  - | b^VAq  + 


4 T 4 T 2 T 9 

+ -,b9  VAq0  + -,q.  B,VAq  + ± q1  A q 

T ~ T T ' 

o o o 


Eq.  (60)  is  used  to  eliminate  q , giving: 
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T T„  Tn  T,  . . 1.  T, 


o = - b3e  + [b^  - b/w^]  qx  + fq1‘[B1  - B^b/]  ^ + 

t 

- | b2TV  J Ex(t-x)  { [VTB2T  - AVTB2T]q1  + V^efdr  + 

O 

[/  Ex(t-t)  ®dx  y*Ex( t-T){vTB2T  - AVTB2T|qldTj 


4 T 

+ ^vi 

o L'o 


+ 


{B2V  - B2VA}Ex(t-T)dT 


T T 

- av  e2  Hidr 


(70) 


2x(t-T){vTB2T  + 


r JL° 

fl 

2 T 

+ T *2  V 
o 

/ Ex(t-x)6dx 

J 

JO 

T 

vb2 


Since  eq.  (70)  must  hold  for  arbitrarily  small  values  of  6 and  q^  , 
the  sum  of  the  first-order  terms  must  itself  approach  zero. 

Thus : 


V + fhTv f 


Ex(t-T)VTb26(r)dT  = [b^  - b2TWTB2T]q1  + 


♦f 


: b2TV  J Ex(t-T)[AVTB2T  - VTB2T]  q1(r)dx 


(71) 


Comparison  of  eq.  (71)  with  eqs.  (65b&c)  shows  that  eq.  (71)  may  be 
expressed  as 


d_ 

dt 


r 

I m(t-x)0(i)dT 


o 


(t-x)  q^(x)dx 


o 


(72) 


A similar  result  is  derived  by  different  means  by  Christensen  [15]. 
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Examination  of  eqs.  (64)  gives  physical  insight  into  the  nature  of 
the  functions  m and  $ . $ , in  eq.  (64a),  is  the  viscoelastic  material 
function  whose  elastic  analogue  is  the  product  of  the  volumetric  coef- 
ficient of  thermal  expansion  and  the  isothermal  bulk  modulus.  In  eq. 
(64b),  m is  the  creep-type  function  corresponding  to  the  heat  capacity 
at  constant  volume  in  the  elastic  case.  From  the  second  law  of  thermo- 
dynamics, we  see  that  this  function  must  be  non-decreasing.  $ is  harder 
to  classify.  The  contribution  of  the  thermal  expansion  will  be  a 
creep- type  function,  while  the  bulk  modulus  will  be  a relaxation-type 
function.  Their  product  will  be  a hybrid  of  sorts,  probably  increasing 
for  some  short  time  and  then  slowly  decreasing  again. 

While  the  heat  capacity,  thermal  expansion,  and  bulk  modulus  have 
all  been  shown  experimentally  to  exhibit  time-dependent  behavior,  such 
effects  are  usually  small  and  the  functions  can  be  sufficiently  well 
approximated  by  step  functions.  Constant  values  for  , a , and  K 

are  reported  in  most  handbooks.  Taking  these  constants  to  be  appropriate 
asymptotic  values,  and  noting  that  isotropy  requires  that  <)>  take  the 
form 


♦ - - oK(1,1,1,0,0,0)T 


(73) 


(where  the  strain  vector,  , is  arranged  in  the  usual  fashion;  the 
three  shear  components  following  the  three  normal  components),  eq.  (72) 
can  be  simplified  to: 


PoCV9 


- oKT  €,  , 
o kk 


(74) 


A 
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Equation  (74)  specifies  the  thermodynamic  relationship  between 
volumetric  strain  and  the  induced  thermal  change.  It  should  be  noted 
that  a very  important  phenomenon,  that  of  heating  due  to  shear  dissi- 
pation, is  not  accounted  for  in  this  formula.  The  contribution  of 
shear  deformation  to  thermal  change  has  been  lost  in  the  lineariza- 
tion which  entered  the  derivation  in  eqs.  (54)  and  (71). 

Heating  due  to  shear  deformation  must  be  independent  of  the 

sense  of  shear  when  the  deformation  occurs  from  an  unstrained  state. 

• • 

Consequently,  the  dependence  of  9 on  e must  be  an  even  function 
in  each  of  the  shear  deformation  components  and  cannot  be  included  in 
a linear  approximation. 

To  include  this  phenomenon  in  the  present  report  would  not  only 
require  retention  of  the  quadratic  terms  currently  present  in  eq.  (70), 
but  would  also  necessitate  the  introduction  of  cubic  terms  in  the 
expansion  of  the  Gibbs'  free  enthalpy  (eq.  (53))  so  that  quadratic 
terms  would  appear  in  the  expression  for  entropy  (eq.  (54c)). 

The  ensuing  analysis  becomes  quite  involved  and  is  not  included 
in  this  report.  Derivation  of  a phenomenological  expression  for  heating 
due  to  shear  dissipation  remains  a subject  for  future  investigation. 


CHAPTER  FOUR 


COMPARISON  AND  SUMMARY 

Although  the  derivations  necessary  to  find  the  temperature  change 
due  to  deformation  differ  in  a great  many  details  between  the  infini- 
tesimal elastic  and  viscoelastic  cases,  the  results  (eqs.  (15)  and  (74)) 
are  essentially  identical.  It  remains  to  evaluate  the  material  proper- 
ties of  significance,  and  to  determine  from  them  the  extent  to  which 
thermomechanical  coupling  effects  need  to  be  considered. 

Since  eq.  (74)  will  be  applied  to  small  deformations  of  metals, 
rubbers,  and  plastics  alike,  data  are  presented  in  Table  Three  for  all 
three  classes  of  material.  Calculations  nave  been  made  for  a temperature 
of  25°C.  ^£jcjc  t^ie  volumetric  strain  necessary  to  induce  a tempera- 
ture change  of  one  degree  centrigrade.  P is  the  hydrostatic  pressure 
necessary  to  cause  that  volumetric  strain.  is  the  longitudinal 

strain  associated  with  Ae^  in  a uniaxial  tension  test.  Units  used 

-5  3 4 3 

in  Table  Three  are:  ct  , 10  /°C  ; p , gram/cm  ; E , K , 10  cal/cm  ; 

Cv  , cal/ gram/ °C  ; P , 10^  psi . 

The  volumetric  strain  associated  with  a one  degree  change  in  tem- 
perature is  seen  to  vary  from  about  0.1%  to  0.7%  , with  the  polymeric 
materials  exhibiting  much  more  diversity  than  the  metals.  The  compres- 
sive stress  associated  with  this  volume  change  varies  from  500  to  4000 
psi  among  the  polymers,  and  from  700  to  5000  psi  among  the  metals. 

Also  calculated  is  the  strain  that  is  necessary  in  uniaxial  tension 
to  cause  a temperature  drop  of  one  degree  centigrade.  The  values  for 
metals  are  fairly  consistant,  ranging  from  0.44%  to  0.60%,  with  the 
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TABLE  3. 

TEMPERATURE 

CHANGE  VS. 

VOLUMETRIC  1 

STRAIN9 

MATERIAL 

a 

K 

E 

po 

CV 

aKT 

P c 

o V 

A£kk 

P 

Elin 

" 

— — 

" 

' 

METALS 

Aluminum 

7.5 

1.80 

1.7 

2.70 

.204 

730. 

.0014 

1.5 

.0044 

Copper 

5.0 

3.29 

3.2 

8.96 

.0893 

610. 

.0016 

3.2 

.0050 

Iron 

3.1 

3.34 

4.4 

7.39 

.107 

390. 

.0026 

5.3 

.0060 

Lead 

8.7 

1.10 

0.40 

11.35 

.0281 

890. 

.0011 

0.73 

.0090 

Nickel 

3.9 

4.49 

5.1 

8.90 

.105 

560. 

.0018 

4.9 

.0048 

Silver 

5.7 

2.48 

1.9 

10.5 

.0540 

740. 

.0013 

2.0 

.0052 

Tin 

6.0 

1.39 

1.2 

7.31 

.0523 

650. 

.0015 

1.3 

.0054 

Zinc 

10.5 

1.66 

2.5 

7.0 

.0851 

870. 

.0011 

1.1 

.0022 

Steel 

3.5 

3.6 

4.5 

7.9 

.118 

400. 

.0025 

5.4 

.0059 

PUBBERS 

Polyisoprene 

67. 

.047 

.913 

.386 

270. 

.0038 

0.11 

(Hevea) 

Polyisoprene 

19. 

.10 

.072 

1.17 

.322 

150. 

.0067 

0.41 

.028 

(Ebonite) 

P butadiene- 

53. 

.060 

.11 

1.15 

.314 

260. 

.0038 

0.14 

.0063 

co-styrene 

Neoprene 

60. 

.050 

1.23 

.476 

150. 

.0066 

0.20 

PLASTICS 

Polystyrene 

23. 

.072 

.079 

1.05 

.28 

170. 

.0060 

0.26 

.016 

Polyvinyl- 

24. 

.098 

1.41 

.24 

210. 

.0048 

0.29 

chloride 

P 'methyl  Meth- 

■ 27. 

.096 

.099 

1.17 

.31 

210. 

.0047 

0.27 

.014 

acrylate 

Low  density 

.079 

.024 

.91 

.25 

780. 

.0013 

0.062 

.013 

P'ethylene 
High  density 

(75.) 

.11 

.96 

.21 

1200. 

.00082  0.055 

P'ethylene 

Nylon  6 

62. 

.079 

1.1 

.26 

510. 

.0020 

0.096 

9 Data  from  [5),  [9],  [44]  for  metals;  [7],  sec.  6.,  for  rubbers;  and 
[8],  tables  4.10,  5.3,  and  9.3  for  polymers. 
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exceptions  of  lead  (0.90%)  and  zinc  (0.22%).  These  two  metals  differ 

primarily  in  their  Poisson's  ratios.  It  should  be  noted  that  e,. 

lin 

exceeds  the  elastic  limit  for  most,  if  not  all,  the  metals  presented. 

This  does  not  imply  that  such  strains  are  obtainable  elastically;  it 
is  merely  a result  of  the  choice  of  1°  C.  as  a benchmark.  Since  the 
strain-temperature  effect  is  linear,  these  values  can  be  used  in  calcu- 
lating the  thermal  effect  of  lesser  strains. 

Calculations  were  also  made  for  those  polymers  for  which  values  of 
Young's  modulus  were  found.  These  materials  show  more  variation; 

ranging  from  0.6%  to  almost  3%  . The  discussion  concerning  e^n 
for  metals  is  applicable  to  the  rubbers  and  plastics  as  well. 

The  thermomechanical  coupling  effect  has  been  shown  to  be  insignif- 
icant in  the  infinitesimal  transient  deformation  of  solids  in  almost 
all  cases.  It  should  be  reiterated  that  eqs.  (15),  (26),  and  (74)  can 
be  used  to  accurately  predict  these  small  thermal  changes,  allowing  one 
to  decide,  in  light  of  the  known  thermal  dependence  of  the  mechanical 
properties  and  the  accuracy  desired,  if  thermomechanical  coupling  effects 
need  be  considered.  It  has  been  shown  in  earlier  chapters  that  the 
equivalent  heat  argument  is  wholly  unsatisfactory  for  such  purposes. 

The  insignificance  of  the  thermomechanical  coupling  effect  has  also 
been  shown  for  the  finite  elastic  deformation  of  rubbers.  Finite,  non- 
linear deformations  of  metals  and  viscoelastic  materials  are  not  con- 
sidered in  this  work. 
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